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ABSTRACT 


A  two-dimensional  numerical  solution  of  the  thermal  Reynolds', 
elasticity,  and  energy  equations  resulted  in  predictions  of  the  flow 
velocities  and  temperatures  within  the  lubricant  film  between  a  crown 
roller  and  a  rotating  shaft.  The  results  showed  a  good  agreement  with 
previously  measured  experimental  data  on  the  surface  temperatures  in  an 
EHD  contact  and  on  the  resulting  frictional  torques  for  both  rolling  and 
sliding.  A  pressure  and  temperature  dependent  viscosity  was  assumed  for 
the  synthetic  paraffinic  hydrocarbon  used  in  the  experiments.  This  model 
allowed  an  evaluation  of  the  relative  importance  of  the  generation, 
viscous  dissipation,  and  compression  terms  in  the  energy  equation  at 
different  stations  along  the  mid-film  line  of  contact.  Based  on  the 
insight  about  the  relative  importance  of  terms  as  found  in  the  numerical 
simulation,  convenient  expressions  were  developed  for  an  approximate 
method  of  predicting  the  mid-film  temperature  and  sliding  friction 
force.  - 


The  influence  of  thermal  expansion  of  the  contact  between  a  crowned 
surface  in  sliding  contact  with  a  rotating  cylindrical  surface  was 
included  in  a  model  of  thermoelastic  instability  to  predict  the  onset  of 
scuffing  failure.  The  simplified  approach  assumed  a  heat  balance  between 
the  generation  from  friction  and  the  heat  conducted  into  the  cylinders.  — 
In  addition  it  was  assumed  that  elastic  deformation  caused  by  the  thermal 
disturbance  had  the  form  of  a  cosine  curve,  and  that  the  heat  conduction 
into  the  rotating  cylinder  was  one-dimensional.  Results  predicting  the 
scuffing  load  at  different  sliding  speeds  showed  reasonable  agreement 
with  previously  reported  experimental  results.  The  scuffing  load  was 
quite  sensitive  to  the  value  of  the  friction  coefficient. 

An  experimental  apparatus  was  built  whereby  the  onset  of  scuffing 
could  be  observed  for  rolling/sliding  contact  between  a  crowned  cylinder 
and  3  straight  cylinder.  Instrumentation  included  the  means  to  measure 
speed,  load,  surface  temperature  near  the  contact  spot,  and  frictional 
drag.  The  onset  of  scuffing  was  sensed  by  increased  vibration  or 
temperature.  Experiments  were  performed  using  three  contact  lubricants 
of  different  rheology,  and  two  steels  with  different  thermal  and  material 
properties.  Results  correlate  with  the  thermoelastic  instability  model 
in  the  regime  of  high  calculated  specific  film  thickness  foi  the  two 
steels  tested. 
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INTRODUCTION 


The  objective  of  the  proposed  program  was  to  analyze  the  thermal 
effects  in  a  rolling/sliding  EHD  contact  to  determine  if  they  represent 
an  important  parameter  leading  to  scuffing  failure.  A  theoretical  model 
determined  the  change  in  shape  due  to  thermal  expansion  of  the 

interacting  surfaces  as  a  result  of  non-uniform  temperatures.  The 
operating  conditions  and  lubricant  properties  were  studied  to  determine 
what  situations  are  unstable;  that  is,  which  lead  to  more  distortion 
and  still  higher  heating. 

Scuffing  is  a  common  mode  of  failure  of  lubricated  contacts  in 
gears,  cam-followers  seals,  and  other  mechanisms  which  have  a  significant 
proportion  of  sliding  in  a  load  cycle.  These  mechanisms  may  be  well 
supplied  with  lubricant,  but  the  film  thickness  will  be  relatively  small 
and  the  heat  generation  will  be  high  because  of  the  sliding.  Scuffing 
failure  occurs  at  small  regions  of  the  apparent  contact  area  and 
typically  consists  of  material  transfer  from  the  stationary  member  to  the 
moving  part.  The  transition  to  scuffing  failure  takes  place  after  an 
initial  period  of  operation  during  which  little  wear  takes  place  and  the 
friction  coefficient  is  low.  Failure  is  exhibited  by  a  rapid  increase  in 
friction  and  high  surface  temperature,  plus  increased  vertical  vibration 
due  to  the  transfer  of  material. 

One  explanation  for  such  a  failure  relates  to  the  thermal  conditions 
in  the  contact  and  their  effect  on  the  lubrication.  The  transition 
period  may  be  a  period  of  heating  and  distorting  the  contact  such  that 
high  pressure  and  temperatures  reduce  the  film  thickness  or  high  tempera¬ 
tures  reduce  or  eliminate  the  effectiveness  of  the  boundary  lubricants 

The  research  program  explored  the  possibility  that  thermal  effects 
in  an  EHD  contact  can  distort  the  surfaces  to  produce  unstable  heating; 
that  is,  local  heating  will  produce  a  change  in  surface  shape  that  will 
lead  to  still  higher  heating  in  that  region.  The  research  involved 
critical  analysis  and  experiments  to  determine  if  such  an  explanation 
could  be  applied  to  scuffing  failures.  The  results  which  demonstrated 
the  validity  of  this  hypothesis  show  that  the  analysis  can  serve  as  a 
basis  for  material  property  and  operating  conditions  selection  to  avoid 
the  occurrence  of  this  failure  mode. 

The  research  is  organized  in  a  series  of  steps  or  tasks  that  will 
ultimately  determine  the  validity  of  the  hypothesis  that  a  thermoelastic 
instability  can  give  rise  to  scuffing.  A  thermal  analysis  of  EHD 
lubrication  in  rolling  and  sliding  contacts  has  been  carried  out  and  is 
described  in  this  report,  as  is  an  elementary  analysis  of  thermoelastic 
instability  for  a  simplified  case.  These  two  steps  have  been  followed  by 
an  experimental  evaluation  and  a  more  detailed  thermoelastic  instability 
analysis  using  the  thermal  EHD  model.  Results  of  the  work  done  are 
presented  in  the  sections  that  follow. 
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SUMMARY 

The  efforts  over  this  reporting  period  have  concentrated  on  three 
aspects.  The  first  was  the  continued  development  of  a  two-dimensional 
numerical  solution  to  the  thermal  problem  of  elastohydrodynamic  lubrica¬ 
tion  of  rolling/sliding  contacts.  The  second  has  been  the  simplified 
formulation  of  a  thermoelastic  instability  theory  applied  to  the  case  of 
thermal  expansion  effects  of  a  stationary  crowned  surface  in  sliding 
contact  with  a  rotating  cylindrical  surface.  The  third  has  been  the 
design  and  testing  of  an  apparatus  to  experimentally  investigate  the 
onset  of  scuffing  and  its  relation  to  load,  speed,  and  local  temperature 
The  overall  goal  is  to  understand  the  role  that  thermal  effects  play  in 
elastohydrodynamic  lubrication  and  eventual  failure  of  the  contact. 


Two-Dimensional  Analysis  of  EHD  Lubrication 

A  two-dimensional  numerical  solution  of  the  thermal  Reynolds' 
elasticity,  and  energy  equations  resulted  in  predictions  of  the  flow 
velocities  and  temperatures  within  the  lubricant  film  between  a  crown 
roller  and  a  rotating  shaft.  The  results  showed  a  good  agreement  with 
previously  measured  experimental  data  on  the  surface  temperatures  in  an 
EHD  contact  and  on  the  resulting  frictional  torques  for  both  rolling  and 
sliding.  A  pressure  and  temperature  dependent  viscosity  was  assumed  for 
the  synthetic  paraffinic  hydrocarbon  used  in  the  experiments.  This  model 
allowed  an  evaluation  of  the  relative  importance  of  the  generation, 
viscous  dissipation,  and  compression  terms  in  the  energy  equation  at 
different  stations  along  the  mid-film  line  of  contact.  Based  on  the 
insight  about  the  relative  importance  of  terms  as  found  in  the  numerical 
simulation,  convenient  expressions  were  developed  for  an  approximate 
method  of  predicting  the  mid-film  temperature  and  sliding  friction 
force . 


The  following  conclusions  were  drawn: 

1.  For  pure  rolling  conditions  the  compression  heating  is  the 
major  contributor  to  the  increase  of  temperature  in  the  contact 
zone . 

2.  There  is  significant  increase  in  film  temperature  as  the  slip 
is  increased. 

3.  The  contribution  of  convection  in  the  y~direction  is  signifi¬ 
cant  in  the  inlet  zone. 

4.  For  slip  conditions  the  major  sources  of  heat  generation  and 
dissipation  are  the  shear  and  conduction. 

5.  For  slip  conditions  as  the  load  increases  the  contribution  of 
compression  and  convection  in  the  x-direction  increase. 

6.  For  slip  conditions  the  maximum  temperature  occurs  at  the 
mid-film  at  a  point  where  the  pressure  is  at  its  greatest  and 
dP/dx  =  0. 
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Approximate  Analysis  for  Estimating  Mid-Film  Temperature 
and  Sliding  Friction  Force 

An  approximate  solution  to  the  problem  of  thermal  elastohydrodynamic 
lubrication  of  sliding  contacts  is  obtained.  An  equation  for  the 
mid-film  temperature  which  is  a  function  of  viscosity,  fluid  conduc¬ 
tivity,  sliding  velocity,  and  temperature  viscosity  exponent  is  obtained. 
The  shear  stress  on  the  disk  surface  over  the  contact  region  is  also 
obtained.  The  following  conclusions  were  drawn: 

1.  The  maximum  temperature  rise  occurs  at  the  middle  of  the  film 
at  the  point  where  dP/dx  =  0. 

2.  The'-e  is  significant  increase  in  the  film  temperature  as  the 
slip  is  increased. 

3.  For  slip  conditions  the  major  sources  of  heat  generation  and 
dissipation  are  the  shear  and  conduction. 

The  mid-film  temperature  increses  linearly  as  the  load  is 
increased. 

5.  For  a  given  load,  the  traction  force  initially  rises  with 
sliding  speed,  reaches  a  maximum,  then  drops  with  further 
increases  in  speed. 

6.  The  results  for  the  mid-film  temperature  and  sliding  traction 
force  obtained  from  FDA  [4]  and  the  approximate  solution 
indicate  excellent  correlation  exists  between  the  results. 


Thermoelastic  Stability  Analysis 

Thermoelastic  instability  in  an  elliptic  contact  between  a  straight 
and  crowned  cylinder  is  suggested  as  a  potential  initiator  of  scuffing 
failure.  An  approximate  equation  is  presented  as  a  guide  to  delineate 
the  region  of  instability.  The  agreement  of  the  calculated  results  with 
the  available  experimental  data  is  satisfactory. 

Since  this  is  a  simplified  approach,  some  assumptions  are  made,  such 
as  the  elastic  deformation  caused  by  the  thermal  disturbance  being  in  the 
form  of  cosine  curve,  and  the  conduction  at  the  surface  of  the 
rotating  cylinder  is  one  dimensional.  Further  theoretical  analysis  is 
being  performed  to  examine  the  thermoelastic  instability  theory. 


Experimental  Tests  of  Scuffing 

An  experimental  apparatus  has  been  built  whereby  the  onset  of 
scuffing  may  be  observed.  Instrumentation  includes  the  measurement  of 
speed,  load,  surface  temperature  near  the  contact  spot,  and  frictional 
drag.  The  onset  of  scuffing  can  be  sensed  by  increased  vibration  or 
temperature.  The  experimental  results  can  be  summarized  as  follows: 


1.  The  theoretical  scuffing  model  is  viable  in  the  regime  of  high 
calculated  specific  film  thickness  (\  >  4)  for  steels  having 
different  thermal  properties. 

2.  The  model  fails  to  adequately  predict  scuffing  in  the  regime  of  low 
calculated  values  of  specific  film  thickness  (X  <  1.3). 

3.  Although  the  theoretical  model  predicts  a  response  to  variation  in 
thermal  and  material  properties,  the  experimental  results  did  not 
answer  the  question  as  to  whether  the  difference  exists  or  not. 

4.  A  time  dependency  to  failure  was  observed  in  the  regime  of  high 

calculated  specific  film  thickness.  The  trend  has  not  been 

quantified  but  is  found  to  agree  with  the  observations  in  other 
research. 

5.  Some  indications  of  moving  hot  spots  were  observed. 
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DETAILS  OF  TECHNICAL  PROGRAM 

During  this  program  three  areas  of  research  have  been  pursued.  The 
first  was  the  continued  development  of  a  two-dimensional  numerical 
solution  to  the  thermal  problem  of  elastohydrcuynaraic  lubrication  of 
rolling/sliding  contacts.  The  second  has  been  the  simplified  formulation 
of  a  therraoelastic  instability  theory  applied  to  the  case  of  thermal 
expansion  effects  of  a  stationary  crowned  surface  in  sliding  contact  with 
a  rotating  cylindrical  surface.  The  third  has  been  the  design  and 
testing  of  an  apparatus  to  experimentally  investigate  the  onset  of 
scuffing  and  its  relation  to  load,  speed  and  local  temperature. 


Thermal  Analysis  of  Elas tohydrodynamic  Lubrication  in  Rolling/Sliding 

Contacts 

A  two  dimensional  numerical  solution  to  the  problem  of  thermal 
elastohydrodynamic  lubrication  of  rolliug/sliding  contacts  was  obtained 
using  a  finite  difference  formulation.  The  technique  involves  the 
simultaneous  solution  of  the  thermal  Reynold's  equation,  and  the  two 

dimensional  energy  equation.  A  pressure  and  temperature  dependent 

viscosity  for  a  synthetic  paraffinic  hydrocarbon  (XRM-109F)  was 

considered  in  the  solution  of  the  Reynolds'  and  energy  equations.  The 

experimental  pressure  and  surface  temperature  measurements  obtained  by 
Dow  and  Kannel  [1]  were  used  in  evaluating  the  results  of  the  numerical 
analysis  for  the  cases  of  pure  rolling  and  slip  conditions. 

Analytical  Development 

Solution  of  the  two  dimensional  thermal  elastohydrodynamic  lubrica¬ 
tion  problem  requires  the  simultaneous  solution  of  the  thermal  Reynold' 
equation,  the  energy  equation  with  the  fluid  properties  as  functions  of 
pressure  and  temperature,  and  the  elasticity  equation. 

The  thermal  Reynolds'  equation  is  obtained  from  the  principles  of 
conservation  of  mass  and  momentum  applied  to  an  element  of  the  lubricant. 
The  following  assumptions  have  been  made  in  the  development  of  the 
governing  equations. 

1.  The  viscosity  of  the  lubricant  is  only  a  function  of  pressure  and 
temperature  (not  dependent  on  the  shear  rate). 

2.  The  body  and  the  inertia  forces  are  negligible  in  comparison  to  the 
pressure  and  viscous  forces. 

3.  The  radii  of  curvature  of  rolling  elements  are  much  larger  than  the 
lubricant  film  thickness.  Therefore,  all  effects  due  to  curvature  of 
the  fluid  film  are  neglected. 

4.  The  pressure  is  constant  across  the  thickness  of  the  film. 

5.  Thermal  conductivity,  specific  heat,  and  thermal  expansibility  do  not 
vary  with  pressure  or  temperature. 
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Using  the  foregoing  assumptions,  the  thermal  Reynolds'  equation,  can 
be  written  as: 

7T  (G2  d7'  =  U1  77'  (G3)  +  ~  u~  (u2  “  Ul)  (1) 

e0 


where 


h  y 

r  r  r  ^  dy 
0  0  u 


c9  =  —  G  -  fhf  fy  dy 
2  1  »  0  ” 


G  =  r  dy 
0 


h 

j  =  f  ly 

eo  o  u 


u  =  r  ^iy 
e !  0  ;J 


The  temperature  field  within  the  lubricant  film  is  obtained  from  the 
solution  of  the  energy  equation.  Neglecting  the  heat  conduction  in  the 
direction  of  flow,  the  energy  equation  may  be  written  as: 

0fcf  ru  4^  +  ^  =  u  +  «Tu  ~  +  k  ~  (7) 
f  f  3x  3y  f  2 


where 


y  u_  -  u,  u  dP/dx  y, 

dP  ry  ,  ,  2  1  e  r  dy  . 

dx  n  u  u  u  _  u  1 

0  e0  e0  _  0 


Ru  =  J_  dP  f 

3y  u  dx  ^  u 


1±)  +  °2  ~  U1 
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The  energy  equation 
(temperature  on  the 
surface  temperature 
for  the  calculation 


(7)  requires  two  boundary  conditions  for  its  solution 
surface  of  each  rolling  element).  The  experimental 
measurements  [1J  were  used  as  the  boundary  conditions 
of  the  temperature  field. 


The  pressure  distributions  from  the  experimental  measurements  [1]  at 
the  center  of  the  the  contact  width  closely  resemble  the  elliptical 
Hertzian  pressure  distribution.  Therefore,  the  inlet  zone  film  shape  is 
assumed  to  be  that  of  the  Hertz  equation  of  dry  line  contact. 


- 

1/2 

2 

1 

-In  r-  + 
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The  Hertz  equation  (10)  is  used  to  calculate  the  film  shape  outside  the 
contact  zone.  The  film  thickness  inside  the  contact  zone  is  calculated 
from  the  thermal  Reynolds'  equation  (1). 

The  viscosity /pressure  relationship  used  in  this  analysis  was 
proposed  by  Roelands  [2]; 

z 

log  u  =  (log  u  +  1.2)  ( 1  +  2 QOQ- ^  -  1*2  (11) 

An  exponential  change  of  viscosity  with  temperature  was  incorporated  and 
the  complete  viscosity  relationship  can  be  expressed  as: 


u  = 


10^log  U0  +  1  +  2000 

*(T  -  T  ) 
e  s 


-) 


1.2] 


(12) 


The  lubricant  under  investigation  was  a  synthetic  parafinic  hydrocarbon. 


The  velocity  profile  in  the  x-direction  obtained  from  the  equili¬ 
brium  of  an  infinitesimal  element  is  given  by  equation  (8).  The 
velocity  profile  in  the  y-direction  (across  the  film)  could  be  obtained 
from  the  continuity  equation: 


Tu  R  v 
3^  +  Ty 


0 


(13) 


This  equation  is  a  first  order  partial  differential  equation  needing  only 
one  boundary  equation  for  its  solution.  However,  the  geometry  of  an  EHD 
regime  illustrated  in  Figure  1  provides  two  boundary  conditions.  If  the 
continuity  equation  is  differentiated  with  respect  to  y  the  following  is 
obtained. 


figure  i.  Velocity  distribution  in  heavily  loaded  contact 
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This  equation  is  a  second  order  partial  differential  equation  needing  two 
boundary  conditions,  which  are  readily  available  from  the  surface 
velocities  or  the  disks.  Equation  (14)  is  used  to  evaluate  velocity  in 
the  y-direction  using  finite  difference  formulation.  Figure  1  illus¬ 
trates  the  variation  of  velocity  in  the  rolling  direction;  backflow 
occurs  in  the  inlet  zone  and  becomes  a  positive  linear  flow  at  the  center 
of  the  contact  where  the  pressure  gradient  is  zero. 

Computational  Scheme 

In  order  to  solve  the  thermal  elastohydrodynamic  lubrication  prob¬ 
lem,  the  isothermal  case  is  solved  first.  The  results  obtained  from  the 
isothermal  case  are  then  used  as  initial  solution  for  the  thermal  case. 
The  isothermal  equations  can  easily  be  obtained  from  the  thermal 
equations  by  allowing  the  viscosity  of  the  lubricant  to  be  a  function  of 
pressure  and  independent  of  temperature. 

The  differential  equations  (1)  and  (7)  were  replaced  by  their 
equivalent  finite  difference  formulation.  An  iterative  solution  of  the 
momentum  and  energy  equations  is  used  to  determine  the  shear  stress,  the 
velocities  in  the  x  and  y-directions ,  and  the  temperature  profile 
throughout  the  inlet  and  contact  zone  for  prescribed  boundary  conditions. 
The  results  of  the  isothermal  case  are  then  used  in  conjunction  with  the 
energy  equation  to  calculate  the  initial  temperature  field  within  the 
lubricant  film.  The  effects  of  this  temperature  distribution  are 
introduced  on  the  viscosity,  film  shape,  and  fluid  flow  in  the  x  and 
y-directions.  These  values  are  then  used  to  calculate  a  new  temperature 
field.  This  iterative  procedure  is  continued  until  the  temperature  field 
reaches  steady  state.  The  rheological  properties  of  the  lubricant  are 
allowed  to  vary  as  a  function  of  local  pressure  and  temperature. 

One  of  the  main  difficulties  in  solving  fluid  dynamic  problems 
involving  the  Navier  Stokes  and  energy  equations,  is  the  appearance  of 
nonlinear  substantial  derivatives.  These  derivatives  require  a  knowledge 
of  the  flow  direction  at  each  node  before  they  can  be  evaluated.  A 
formalized  method  of  handling  such  problems  is  to  utilize  the  "upwind" 
differencing  scheme  [3].  The  temperature  at  each  surface  node  is  known 
from  measured  surface  temperatures.  It  is  also  assumed  to  be  ambient  at 
the  entry  nodes  where  the  x-component  of  velocity,  u  is  positive.  If  the 
flow  field  predicts  backflow  at  some  entry  nodes,  the  temperature  at 
these  nodes  are  dependent  values  and  must  be  found  during  the  numerical 
solution.  The  major  difference  between  the  "upwind"  differencing  scheme 
and  other  solution  schemes  lies  in  the  use  of  upwind  values  to  define  a 
convect  flux.  The  local  velocities  in  the  x  and  y-directions  are  defined 
from  equations  (8)  and  (14).  If  the  flow  is  positive  the  downwind  value 
for  the  temperature  at  node  (i,j)  is  assumed;  and  if  the  flow  is  negative 
the  upwind  value  is  assumed  for  that  node.  A  number  of  subroutines  have 
been  constructed  for  the  evaluation  of  the  temperature  field  within  the 
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lubricant  film.  These  subroutines  provide  information  such  a  transient 
viscosity,  viscosity  integrals,  etc.  Values  are  recalculated  using  these 
subroutines  after  each  temperature  field  is  obtained.  The  fluid  inte¬ 
grals  (equations  (2),  (3),  (5),  and  (6))  were  evaluated  using  the 
trapezoidal  rule. 

The  flow  chart  of  Figure  2  illustrates  the  computational  procedure 
that  was  used  for  calculating  the  temperature  field  within  the  lubricant 
film.  The  isothermal  viscosity  of  the  lubricant  (equation  (11)),  film 
shape  (equation  10),  and  velocity  in  the  x  and  y-direction  (equations  (8) 
and  (14))  are  calculated  based  on  the  measured  pressure  and  surface 
temperature  data.  These  values  are  used  to  arrive  at  the  initial 
temperature  field.  The  calculated  temperature  field  is  then  used  to 
introduce  the  influence  of  temperature  on  the  viscosity  and  a  new  flow 
field  is  calculated.  The  iteration  procedure  is  continued  until  the 
lubricant  temperature  obtained  from  the  energy  equation  (7)  does  not 
change  from  one  iteration  to  the  next.  One  aspect  of  the  solution 
technique  is  the  limitation  on  the  temperature  change  per  iteration. 
This  procedure  was  designed  to  speed  convergence  and  involves  limiting 
the  temperature  change  which  can  occur  for  a  flow  field  iteration.  The 
final  temperature  distribution  is  the  same  with  or  without  using  a 
limited  temperature  change,  but  convergence  is  reached  considerably 
faster  when  the  limiting  temperature  is  used.  Using  the  limitation  that 
the  temperature  change  will  be  less  than  or  equal  to  10%  of  the  current 
temperature,  the  energy  equation  can  be  solved  in  15-20  iterations. 

Results  and  Discussion  of  Numerical  Solution 


Figures  3  and  4  illustrate  the  experimental  pressure  and  surface 
temperature  measurements  obtained  by  Dow  and  Kannel  [1]  using  a  rolling 
disk  machine  with  a  synthetic  paraffinic  oil.  Figures  3a  and  3b  depict 
pressure  and  surface  temperature  for  pure  rolling  and  6.2%  slip  condition 
for  a  disk  load  of  970N.  The  dimensions  and  the  properties  of  the 
rolling  disk  machine  and  lubricant  are  tabulated  in  Table  1.  Figure  3a 
indicates  that  the  peak  surface  temperature  for  pure  rolling  occurs  in 
the  inlet  zone  near  the  point  of  maximum  pressure  gradient.  For  the  case 
of  6.2%  slip  the  peak  surface  temperature  has  moved  toward  the  center  of 
the  contact  where  the  pressure  is  maximum.  Figure  4  illustrates  the 
meaaured  pressure  and  surface  temperature  at  454  N  of  disk  load  for  pure 
rolling  and  5.3%  slip  condition.  For  the  pure  rolling  (Figure  4a),  peak 
surface  temperature  rise  is  12°  C.  above  the  ambient  and  occurs  in  the 
inlet  zone  near  the  peak  pressure  gradient  as  it  did  of  for  the  970N 
load.  For  the  case  of  5.3%  slip  (Figure  4b),  the  temperature  peak  has 
again  moved  towards  the  center  of  the  contact  to  the  point  of  maximum 
pressure  as  in  the  case  of  6.2%  slip  at  970  N  disk  load.  These 
experimental  measurements  were  used  as  the  boundary  conditions  for  the 
solution  of  thermal  Reynolds’  equation,  and  the  two  dimensional  heat 
equat ion. 

Figures  5  and  6  illustrate  the  calculated  temperature  distribution 
through  the  thickness  of  the  lubricant  film  resulting  from  the  effects  of 
shear,  compression,  conduction,  and  convection  in  the  x  and  y-directions. 
For  pure  rolling  condition  (Figure  5a),  the  peak  lubricant  temperature 
rise  in  the  film  is  36°  C.  above  the  ambient  whereas  the  peak  surface 


Figure  3.  Dimensionless  pressure  and  surface  temperature  variation  through  the  center 
of  contact  (970  N). 


Direction,  Length  (cm) 
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Table  1.  Disk  and  lubriant  properties 

Material  of  both  disks  steel 

Lubricant  synthetic  paraffinic  hydrocarbon  XRM-109F 

Diameters  of  upper  disk,  cm  3.569 

Diameter  of  lower  disk,  cm  3.546 

Crown  radius  of  lower  disk,  cm  20.569 

Rotational  speed  (nominal),  rad/sec  525 

Surface  speed  (nominal),  cm/sec  934.7 

Specific  heat  of  the  disk  3723 

times  the  density,  KJ/(m3  *  C) 

Thermal  conductivity  of  the  27.22 

disk,  W/(m  *  C) 

Modulus  of  elasticity  of  the  207 

disks,  GPa 

Poisson’s  ratio  of  the  disks  .3 

Thermal  conductivity  of  the  240 

lubricant,  W/(m  °  C) 

Specific  heat  of  the  lubricant  1744 

times  the  density,  KJ/(m3  "  C) 

Thermal  viscosity  coefficient,  1/  ”C  .04 

Thermal  expansivity  of  the  lubricant,  1/  "C  .00069 

Exponent  for  Roeland’s  viscosity  model  .45 


Figure  5.  Temperature  distribution  through  the  thickness  o 
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temperature  is  13°  C.  above  the  ambient.  For  the  6.2%  slip  case  Figure 
5b,  the  peak  temperature  in  the  film  is  110°  C.  above  the  ambient  for  a 
surface  peak  of  24°  C.  above  the  ambient.  It  is  important  to  note  that 
the  lubricant  film  temperature  rise  for  the  pure  rolling  condition  mainly 
occurs  in  the  inlet  zone.  Whereas  for  the  slip  condition  the  rise  occurs 
in  the  contact  zone.  The  temperature  rise  for  pure  rolling  condition  is 
moderate,  but  significant  temperature  increase  occurs  for  the  slip 
condition.  Figure  6  depicts  the  temperature  distribution  for  pure 
rolling  and  5.3%  slip  condition  at  454  N  of  disk  load.  Again  in  this 
case,  as  in  the  case  of  970  N  disk  load,  significant  temperature  rise 
occurs  at  the  mid-film  for  the  5.3%  slip  condition.  The  mid-film 
temperature  (Figure  6b)  rises  to  72°  C.  above  ambient  compared  to  the 
surface  temperature  rise  of  16°  C.  above  ambient. 

Figures  7  through  11  illustrate  the  component  of  heat  generation  and 
dissipation  for  pure  rolling  and  6.2%  slip  conditions  at  a  disk  load  of 
970  N.  Close  study  of  these  components  reveals  the  significance  of  each 
and  the  location  of  its  maximum. 

Figure  7  illustrates  the  component  of  heat  generation  due  to  shear. 
For  the  pure  rolling  condition  (Figure  7a)  all  of  the  shear  heating  is 
generated  in  the  inlet  zone  near  the  disk  surface.  The  shear  heating  is 
maximum  at  the  disk  surface  and  gradually  decreases  to  its  minimum  of 
zero  at  the  middle  of  the  film.  The  shear  heating  is  zero  in  the  contact 
zone;  because  the  gradient  of  velocity  is  zero  in  this  region  for  pure 
rolling  conditions.  The  maximum  value  of  shear  heating  occurs  at  the 
surface  of  the  disks  near  the  entrance  to  the  contact  zone.  Figure  7b 
depicts  the  shear  heating  for  6.2%  slip.  For  this  case,  in  contrast  of 
pure  rolling,  most  of  the  shear  heating  occurs  in  the  high  pressure 
region  (contact  zone).  The  maximum  value  is  22  times  greater  than  that 
of  the  pure  rolling  condition.  Significant  shear  heating  also  occurs  in 
the  inlet  zone,  but  compared  to  the  peak  value  it  does  not  appear  large 
on  the  graph. 

Figure  8  illustrates  the  magnitude  and  distribution  of  the  heat 
generated  by  compression  in  the  lubricant  film.  Compressive  heating  is 
proportional  to  the  pressure  gradient,  and  as  a  result  it  will  be  the 
same  across  the  thickness  of  the  film.  The  heating  will  also  be  similiar 
for  pure  rolling  and  slip  conditions;  because  of  their  similiar  pressure 
distributions.  Comparing  the  magnitude  of  the  heat  generated  by  shear 
(Figure  7)  and  compresssion  (Figure  8)  indicates  that  shear  heating  plays 
a  dominant  role  in  the  temperature  rise  across  the  fluid  film  for  slip 
conditions  but  is  similiar  to  compressive  heating  for  pure  rolling 
conditions . 

The  component  of  heat  loss  due  to  conduction  is  presented  in 
Figure  9.  Conduction  heat  loss  is  proportional  to  the  second  derivative 
of  temperature  across  the  film  as  well  as  the  coefficient  of  thermal 
conductivity  of  the  lubricant.  It  is  instructive  to  compare  the  heating 
due  to  shear  (Figure  7)  and  that  due  to  compression  (Figure  8)  with  the 
heat  loss  due  to  conduction  shown  in  Figure  9.  The  peak  value  of 
conduction  for  pure  rolling  is  essentially  all  the  heat  generated  due  to 
compression.  For  this  case  the  shear  heating  is  small.  Figure  9b 
depicts  the  conduction  loss  for  the  6.2%  slip  condition.  As  in  the  case 


b)  6.2%  slip  condition 
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Figure  10.  Heat  loss  due  to  convection 
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of  pure  rolling  condition,  conduction  is  the  dominant  mode  of  heat 
transfer;  however,  for  this  case  the  conduction  heat  loss  in  the  contact 
zone  appears  identical  to  the  shear  heating  gain  in  the  film  (Figure  7b) 
since  the  compressive  heating  is  negligibly  small. 

Convection  heat  loss  in  the  x-direction  (rolling  direction)  is 
proportional  to  the  fluid  density,  specific  heat,  velocity  in  the 
x-direction,  and  the  gradient  of  temperature  in  the  x-direction.  The 
convection  peaks  at  the  mid-film  in  the  center  of  the  contact,  where  the 
velocity  and  the  gradient  of  temperature  are  the  greatest.  Figure  10 
illustrates  the  convection  heat  loss  in  the  x-direction.  Even  though  for 
pure  rolling  (Figure  10a)  the  maximum  convection  is  1/4  of  that  due  to 
conduction;  its  contribution  cannot  be  neglected.  Figure  10b  depicts  the 
convection  in  the  x-direction  for  6.2%  slip.  Its  magnitude  is  much 
larger  than  those  of  pure  rolling  as  the  change  in  scale  indicates. 

Figure  11  illustrates  the  component  of  heat  loss  due  to  convection 
in  the  y-direction  (across  the  film).  This  convection  is  also 
proportional  to  fluid  density,  specific  heat,  velocity  in  the  y- 
direction,  and  the  gradient  of  temperature  in  the  y-direction.  The 
velocity  in  the  y-direction  has  its  maximum  at  the  boundaries  and 
gradually  decreases  to  zero  at  the  mid-film  for  pure  rolling,  and  at  a 
point  close  to  mid-film  for  slip  conditions.  This  convection  term 
reaches  its  maximum  of  (Figure  11a)  at  the  boundary  in  the  inlet  zone 
where  the  velocity  in  the  y-direction  is  the  greatest.  It  is  zero  at  the 
mid-film  and  in  the  contact  zone  where  the  boundaries  are  nearly 
parallel.  Figure  lib  depicts  the  convection  in  the  y-direction  for  6.2% 
slip.  It  has  a  maximum  value  at  the  boundary  as  in  the  case  of  pure 
rolling  and  goes  to  zero  at  a  point  near  the  middle  of  the  film.  Close 
study  of  this  term  reveals  its  contribution  is  only  in  the  inlet  zone, 
and  has  no  effect  on  the  temperature  in  the  contact  zone  and  can  thus  be 
neglected . 


An  Approximate  Method  for  Estimating  Mid-Film  Temperature  and 
Sliding  Friction  Force  in  EHD  Lubrication 

An  approximate  method  for  estimating  mid-film  temperature  and 
sliding  friction  force  in  elastohydrodynamic  (EHD)  lubrication  of  sliding 
contacts  has  been  obtained.  The  calculation  is  an  approximation  based  on 
the  conditions  at  the  point  of  maximum  pressure  where  dP/dx  =  0,  but  the 
results  indicate  excellent  correlation  exists  between  this  approximate 
method  and  results  of  the  more  comprehensive  analysis  presented  in  the 
last  section. 

Basis  for  the  Approximation 


The  results  of  the  last  section  (to  be  published  in  ASME  Journal  of 
Tribology  [4])  can  be  used  to  calculate  the  temperature  in  an  EHD  fluid 
film.  This  solution  takes  into  account  the  conduction,  convection  in 
both  x  and  y-directions ,  shear  and  compression.  It  illustrates  the 
temperature  rise  in  sliding  conditions  and  extensively  investigates  the 
contribution  of  each  component  of  the  heat  generation  and  dissipation. 
The  results  obtained  indicated  that  the  maximum  temperature  rise  occurs 
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at  the  middle  of  the  film.  Therefore,  it  is  of  interest  to  compare  the 
magnitude  of  heat  generation  and  dissipation  at  the  middle  of  the  film. 
Figure  12  represents  these  terms  at  the  middle  of  the  film  for  pure  and 
6.2%  slip  condition  at  a  disk  load  of  9  ton.  For  the  pure  rolling 
condition  (Figure  12a)  it  is  observed  that  the  convection  in  the 
y-direction  and  shear  have  no  influence  on  the  heat  balance  in  the  middle 
of  the  film.  The  major  contributors  are  the  compression  and  conduction. 
Although,  convection  in  the  x-direction  removes  heat  from  the  contact; 
its  influence  is  much  smaller  than  that  of  conduction.  Figure  12b 
depicts  the  generation  and  dissipation  terms  at  the  same  disk  load  at 
6.2%  slip.  In  this  case,  in  contrast  to  pure  conditions,  the  major 
contributor  to  the  temperature  rise  in  the  film  is  shear  heating.  The 
major  loss  term  is  the  conduction.  Though,  the  compression  and  convec¬ 
tion  in  the  x-direction  are  not  negligible,  as  a  first  approximation  only 
the  shear  and  conduction  need  be  considered. 


Figure  13  illustrates  the  heat  generation  and  dissipation  at  454  N 
disk  load  for  pure  rolling  and  5.3%  slip  condition.  This  figure  confirms 
the  conclusion  made  from  the  higher  load  condition  previously  discussed, 
namely,  for  slip  conditions  the  shear  is  the  major  heat  source  and  the 
conduction  is  the  major  heat  sink.  However,  the  contribution  of 
compression  and  convection  in  the  x-direction  are  less  significant  in 
this  load  than  at  the  970  N  load.  Therefore,  the  components  of 
compression  and  convection  in  the  x  and  y-direction  in  the  middle  of  the 
film  could  be  neglected.  Thus,  the  predominant  modes  of  heat  generation 
and  dissipation  are  the  shear  and  conduction.  This  conclusion  is  used  as 
the  basis  for  an  approximate  solution  to  the  problem  of  thermal  EHD 
lubrication. 

Governing  Equations  for  the  Approximate  Analysis 

Assuming  that  the  predominant  terms  in  the  film  are  the  shear 
heating  and  conduction  heat  loss,  the  energy  equation  (1)  can  be  obtained 
from  the  heat  balance  on  an  elemental  control  volume. 
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If  compression  heating  and  convection  of  heat  along  and  across  the  film 
(x  and  y-directions  resnectively )  can  be  neglected,  then  the  energy 
equation  for  a  Newtonian  lubricant  model  can  be  written  as: 
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Rewriting  in  terms  of  shear  stress,  equation  (2)  may  be  written  as: 


(3) 


a)  Pure  rolling  condition  b)  6.2%  slip  condition 
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The  Roelands'  viscosity  model  [6]  including  the  temperature  effect  is 
given  by: 


u 


l0r(log  u 0  +  1.2)0  +  P/ 28400) 2  -  1.21 


(A) 


For  convenience  the  numerator  will  be  represented  by  u(P).  Substituting 
for  viscosity  in  equation  (3)  results  in: 
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This  equation  can  be  easily  integrated  with  respect  to  y  if  the  shear 
stress  (t)  is  constant  with  respect  to  y.  This  condition  exists  at  the 
peak  pressure.  It  was  observed  from  the  analysis  that  the  maximum 
temperature  occurs  at  the  point  of  maximum  pressure,  where  dP/dx  =  0. 
The  basic  relation  for  the  equilibrium  of  an  element  in  an  EHD 
lubrication  contact  is  given  by  equation  (6). 
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Evaluating  this  equation  at  the  point  of  maximum  pressure  results  in: 
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Therefore,  t  is  constant  with  respect  to  y.  The  rest  of  the  analysis 
presented  here  is  thus  only  valid  at  dP/dx  =  0.  Integrating  equation  (5) 
the  following  is  obtained. 
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In  order  to  evaluate  C[,  define  the  position  where 


=  0  at  y  =  yc  (10) 

Where  0<yc<h,  also  at  y  =  yc,  T  =  Tc,  and  u  =  uc.  Hence  the  constant 
is : 
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Substituting  cj  in  equation  (9)  results  in; 
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rearranging  the  above  equation  results  in  the  following. 
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Differentiating  equation  (13)  with  respect  to  T  results  in: 
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Substituting  equation  (17)  into  equation  (16)  and  integrating: 
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The  boundary  conditions  are;  T  =  Tj  at  y  =  0  and  T  =  T2  at  y  =  1 
Tc  at  y  =  yc.  Integrating  equation  (18)  between  y  =  0  and  ; 
following  is  obtained. 
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Substituting  A  and  w  in  equation  (19)  and  assuming  both  surface 
temperatures  to  be  equal  to  Tg ,  the  shear  stress  expression  becomes. 
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Equation  (21)  describes  the  shear  stress  present  in  the  contact.  The 
friction  force  is  obtained  by  multiplying  this  shear  stress  by  the  area 
of  the  contact.  The  length  of  the  contact  area  for  the  approximate 
solution  is  the  length  of  the  high  pressure  regions,  where  the  shear 
stress  is  determined.  Equation  (21)  requires  an  estimate  of  mid-film 
(y-direction)  temperature  therefore,  this  temperature  is  developed  next. 

From  the  Newtonian  lubricant  model  assumption  the  velocity  gradient 
could  be  written  as: 
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Rearranging  the  temperature  gradient  equation  (12) 
obtained. 
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the  following  is 
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Combining  equations  (22)  and  (23)  to  eliminate  dy : 
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Integrating  this  equation  yields: 
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Using  the  boundary  conditions  at  y  =  0,  u  =  uj,  T  =  Tj,  and  at  y  -  h,  u  = 
u2  and  T  =  T2,  results  in: 
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This  equation  allows  Tc  to  be  calculated  in  terras  of  the  boundary 
conditions  U2  and  u^.  Letting  =  ^2  =  equation  (26).  The 
following  is  obtained. 
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Equation  (27)  predicts  the  temperature  at  the  mid-film  at  the  point  of 
maximum  pressure  and  dP/dx  =  0.  Figure  14  compares  the  mid-film 

temperature  above  ambient  at  dP/dx  =  0,  obtained  from  the  finite 

difference  analysis  (FDA)  [4]  and  equation  (27)  as  a  function  of  maximum 
Hertz  pressure  for  different  slip  ratios.  It  illustrates  that  as  the 
pressure  increases  (load  increases),  the  maximum  film  temperature 
increases  linearly.  Figure  14  illustrates  that  the  approximate  tempera¬ 
ture  equation  (27)  predicts  values  of  mid-film  temperature,  which  are  4% 
higher  than  those  obtained  from  FDA  [4].  This  high  prediction  is  mainly 
due  to  the  neglect  of  convection  in  the  x-direction.  It  is  also  apparent 
that  good  correlation  exists  between  the  results  obtained  from  FDA  [4] 
and  the  approximate  solution  of  equation  (27). 

The  mid-film  temperature  equation  (27)  depends  on  the  temperature  at 
the  surface  of  the  disk.  However,  the  surface  temperature  Tg,  is 
typically  unknown.  Table  2  contains  the  mid-film  temperature  obtained 
from  equation  (27)  when  surface  temperature  was  retained  and  when 
neglected.  Table  2  indicates  that  the  contribution  of  Ts  in  equation 
(27)  may  be  neglected  thus  reducing  to: 
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Table  2  also  contains  the  mid-film  temperature  obtained  from  the  FDA  [4], 
The  comparison  indicates  that  equation  (28)  does  not  correlate  well  with 
the  results  from  FDA  [4]  and  equation  (27)  at  low  pressures  and  low  slip 
conditions.  However,  as  the  pressure  is  increased  at  a  particular 
sliding  velocity  the  predictibi lity  of  equation  (28)  greatly  improves. 
The  accuracy  also  increases  as  sliding  velocity  is  increased. 


The  Approximate  Analysis  Compared  with  the  More  Complete  Numerical 
Analysis 


The  variation  of  sliding  friction  force  equation  (21)  with  pressure 
is  illustrated  in  Figure  15.  The  increase  in  pressure,  results  in  an 
increase  in  sliding  force  at  a  constant  sliding  velocity.  This  is  a 
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Figure  14.  Peak  Hertz  pressure  versus  mid-film  temperature 
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Table  2.  Results  obtained  from  equation  (27 ),  equation  (28  ),  and  FDA  (4) . 


Pressure 

GPa 

Equation 

(27) 

3%  slip 

C 

Equation 

(28) 

39o  slip 

C 

FDA  14] 

3  9o  slip 

C 

Equation 

(27) 

6  9o  slip 

C 

Equation 

(28) 

69c  slip 

C 

FDA  [4] 
69c  slip 
C 

.8619 

17.73 

-3.60 

18.22 

40.87 

32.70 

41.94 

1.034 

34.79 

24.42 

35.64 

64.35 

60.72 

64.33 

1.207 

55.32 

50.36 

56.00 

88.32 

86.66 

86.84 

1.379 

77.01 

74.59 

76.36 

111.68 

110.89 

110.23 

1.551 

98.61 

97.41 

98.04 

134.09 

133.70 

131.05 

1.724 

119.63 

119.01 

116.89 

155.51 

155.31 

151.04 

1.896 

139.91 

139.58 

135.90 

175.98 

175.87 

172.87 

1.965 

147.80 

147.54 

144.16 

183.92 

183.84 

179.92 

Pressure 

Equation 

(27) 

Equation 

(28) 

FDA  [4] 

Equation 

(27) 

Equation 

(28) 

FDAL4] 

99o  slip 

9  9o  slip 

9  9a  slip 

159c  slip 

15*70  slip 

159c  slip 

GPa 

C 

C 

C 

C 

C 

C 

.8619 

56.13 

51.33 

56.76 

79.15 

76.88 

79.15 

1.034 

81.43 

79.36 

81.81 

105.87 

104.90 

103.45 

1.207 

106.23 

105.29 

104.19 

131.27 

130.83 

128.22 

1.379 

129.98 

129.52 

126.34 

155.28 

155.06 

150.77 

1.551 

152.56 

152.34 

148.78 

177.99 

177.88 

176.17 

1.724 

174.06 

173.94 

170.92 

199.54 

199.48 

195.30 

1.896 

194.57 

194.51 

189.31 

220.08 

220.05 

217.22 

1.965 

202.52 

202.47 

198.11 

228.04 

228.01 

224.65 
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linear  increase.  The  sliding  friction  force  also  increases  as  the  slip 
ratio  is  increased.  The  error  between  the  approximate  method  obtained 
from  equation  (21)  and  FDA  [4]  seems  to  decrease  as  the  sliding  velocity 
is  increased.  The  results  from  the  FDA  [4]  and  an  approximate  solution 
correlate  very  well  as  illustrated  in  Figure  15.  The  maximum  error 
between  the  FDA  [4]  and  the  approximate  solution  is  14%. 

Figure  16  depicts  the  sliding  friction  force  obtained  from  the  FDA 
[4]  and  equation  (21)  versus  sliding  velocities  at  different  pressure 
(loads).  The  measured  sliding  friction  force  at  pressures  of  1.04  GPa 
and  1.36  GPa  at  approximately  6%  slip  obtained  by  Dow  and  Kannel  [1]  are 
also  depicted.  There  exists  good  agreement  between  the  approximate 
friction  force,  the  measured  friction  force,  and  the  friction  force 
obtained  from  the  FDA  [4].  Equations  (21)  and  (28)  could  be  used  as  a 
first  step  in  studying  the  effects  of  temperature  in  EHD  lubrication  of 
sliding  contacts. 


Thermoelastic  Instability  in  Elliptic  Contact  Between  Cylinders  Under 

Sliding  Contact 

The  scuffing  failure  of  a  sliding  contact  is  a  complicated  process 
influenced  by  surface  physics,  chemistry,  and  mechanical  science.  No 
comprehensive  theory  has  yet  been  presented;  however,  much  has  benn 
written  on  this  mode  of  failure.  The  objective  of  this  (section)  is  to 
present  a  mechanical  hypothesis  of  scuffing  failure  relating  the  heat 
generated  due  to  friction  to  the  resulting  distortion  of  the  mating 
surfaces.  From  the  analysis  presented,  the  conditions  for  the  heating/ 
distortion  process  to  be  unstable  can  be  calculated.  The  analysis 
relates  the  influence  of  load,  speed,  coefficient  of  friction,  and 
material  properties  to  the  occurrance  of  unstable  thermal  conditions 
which  could  lead  to  scuffing  failure.  The  results  from  this  analysis  are 
shown  to  provide  satisfactory  agreement  with  published  experimental  data. 

Conceptual  Framework 

The  frictional  heating  of  two  bodies  in  contact  causes  thermoelastic 
deformation  which  may  concentrate  the  heating  and  greatly  increase  the 
contact  pressure  and  temperature.  This  process  can  lead  to  an  unstable 
feedback  loop  driven  by  frictional  heating  and  thermal  expansion  and 
ameliorated  by  heat  transfer  and  wear.  When  this  process  leads  to  an 
unstable  increase  in  contact  stress,  it  is  known  as  thermoelastic 
instability  (TEI).  Modern  demands  for  increased  performance  of  sliding 
systems  make  the  study  of  thermoelastic  instability  important. 

The  idea  of  a  high  "flash"  temperature  in  a  frictionally  heated 
contact  was  introduced  by  Blok  [5]  and  applied  by  him  to  explain  the 
failure  of  lubricated  Hertzian  contacts  as  encountered  in  gearing.  Holm 
[6],  Archard  [7],  and  others  have  drawn  upon  similar  concepts  to  provide 
an  estimate  of  localized  high  temperatures  in  frictional  sliding. 


Sibley  and  Allen  [8]  carried  out  a  series  of  experiments  on  seal 
materials  and  developed  a  criterion  for  thermal  checking  and  showed 
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Figure  16.  Sliding  velocity  versus  traction  force 
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photographic  evidence  of  systematically  moving  hot  patches  in  the  contact 
zone.  Interest  in  contact  instability  and  patch  formation  was  acceler¬ 
ated  by  the  work  of  Barber  (9-11]  who  has  demonstrated  the  phenomenon 
experimentally  and  has  provided  analysis  which  partially  explain  his 
observations  as  well  as  making  fundamental  contributions  to  the  fields  of 
thermoelasticity.  His  explanations  draw  upon  the  modification  of  asper¬ 
ity  contact  by  frictional  heating  and  wear.  His  initial  interest  was  the 
explanation  of  hot  spot  effects  in  railroad  brakes. 

Dow  and  Burton  [12,13]  addressed  the  2-D  problem  of  a  scraper 
sliding  perpendicular  to  its  edge  on  a  conductive  slab  and  showed  that 
instabilities  could  be  predicted  with  or  without  wear.  Dow  carried  out  a 
numerical  simulation  [14]  which  predicted  the  formation  and  translation 
of  isolated  patches  and  experiments  which  vividly  displayed  patch 
formation  [13].  The  occurrence  of  such  patches  is  close  to  his 
theoretical  critical  sliding  speed  for  dry  contact,  but  does  not 
correctly  predict  the  experimental  results  obtained  for  wick-lubricated 
contact.  Burton,  et.  al  [16]  also  demonstrated  non-uniform  contact 
exists  in  seal-like  configuration. 

Durkee  and  Cheng  [17,18]  studied  the  mechanism  of  scuffing  failure, 
especially  the  analytical  modelling  of  the  rough  elastohydrodynamic 
problem.  They  obtained  valuable  experimental  results  concerning  the  load 
at  failure  and  the  sliding  speed  between  the  contact  of  one  crowned 
roller  and  one  cylinder.  Sovak  and  Cheng  [19]  extended  the  work  to 
examine  the  effect  of  contact  geometry  on  the  nature  of  the  initial 
damage. 

From  the  above  survey  it  can  be  seen  that  great  progress  has  been 
made  over  the  years.  However,  there  is  still  controversy  and  a  number  of 
unanswered  questions  for  research.  The  subtlety  of  the  interactions  and 
the  difficulties  of  experimentation  are  such  that  the  explanation  of  many 
important  effects  is  not  generally  well  known  and  guidelines  are  not  yet 
available  to  the  designer.  This  investigation  is  to  establish  the  theory 
of  thermoelastic  instability  in  the  case  of  elliptic  contact  of  sliding 
cylinders  and  to  present  a  simple  formulation  of  the  problem  to  ascertain 
if  it  can  be  used  as  an  estimate  for  the  onset  of  scuffing  failure. 

Theory  and  Analysis  for  a  Simplified  Case 

The  influence  of  thermal  expansion  on  the  contact  between  crowned 
cylinders  is  investigated.  As  shown  in  Figure  17,  the  upper  cylinder 
with  a  crowned  surface  is  stationary  and  the  lower  one  with  the 
cylindrical  surface  is  rotating. 

The  Thermally  Unstable  Condition  -  Surface  irregularities  or  debris 
may  disturb  the  pressure  distribution  by  raising  the  pressure  at  that 
point  momentarily  and  cause  perturbations  of  frictional  heat  input  and 
temperature  in  the  contact  area.  The  growth  or  decay  of  the  disturbance 
depends  upon  heat  conduction,  thermal  deformation  and  wear.  Wear  may 
have  some  stabilizing  effect,  but  only  heat  conduction  and  thermal 
deformation  will  be  considered  here.  If  the  frictional  heat 
generation  caused  by  the  pressure  perturbation  is  greater  than  the  heat 
conducted  into  the  cylinders,  heat  will  accumulate  within  the  contact, 
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then  the  surfaces  of  the  two  cylinders  will  deform  further,  and  the 
temperature  will  continue  to  rise.  The  heat  balance  due  to  the  pressure 
perturbation  at  the  interface  is: 

Qi+Q2  =  Qf  ( 1 ) 

where  0;  is  the  heat  conducted  into  tne  stationary  cylinder,  Q2  is  the 
heat  conducted  into  the  rotating  cylinder  and  Of  is  the  heat  generated  by 
friction  per  unit  time. 


Heat  Conducted  by  the  Stationary  Cylinder  -  The  heat 
equation  is:  ” 


conduction 


9  2  2 

3  T  ,  3  T  3  T  1  3  T 

22+2k3t^  ^ 

3x  3y  3  z  1 

with  boundary  conditions, 

T  =  0  at  x  =  ±  a/2 

T  ■*  0  as  y  *  “ 

T  =  0  as  z  =  t  b/2 

and  equation  (1)  which  is  the  boundary  condition  at  y  =  0. 

With  reference  to  the  experiments  of  Barber  [10]  and  Sibley  [3],  a 
temperature  disturbance  of  the  following  form  is  selected, 

*t 

T  =  T  F(y)  e  cos  us.x  cos  m.z 
o  1  i. 

The  y-dependence  of  the  temperature  can  be  found  by  substituting 
this  perturbation  into  the  conduction  equation.  The  temperature  of  the 
cylinder  is  then  of  the  form 

y  +  St 

T  =  T  e  cos  u,x  cos  w.z  (3) 

o  12 


where 
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the  conduction  in  the  y-direction  per  unit  area  per  unit  time  is 
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The  heat  input  into  the  surface  at  that  part  of  the  contact  area 
where  the  temperature  is  increased  can  be  found  by  integration  over  this 
area,  i.e. 


0,  -  rr  , 


[y=o 


St 

dxdz  =  K.8  T  e  I. 

1  o  1 
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where 


I  =  f  f  cos  m^x  cos  dxdz 


Heat  Carried  Away  by  the  Moving  Surface  -  The  thermal  problem  of  a 
heat  source  moving  along  a  semi -inf inite  surface  has  been  investigated  by 
Tseng  and  Burton  [20],  Barber  and  Martin-Moran  [21]  and  recently  by  Gecim 
and  Winer  [22].  The  following  derivation  is  a  straightforward  method  to 
get  an  approximate  expression  for  the  heat  input. 


Since  the  surface  speed  is  high  and  the  contact  area  is  small,  the 
heat  conduction  will  be  considered  as  one  dimensional.  The  amount  of 
heat  passing  through  dxdy,  a  small  element  of  the  elliptical  contact 
shown  in  Figure  18,  in  unit  time  [12]  is 

dQ,  =  K  T  - ""T/V  dxdz  ^ 

(^k2t)l/Z 

The  time  of  contact  t  can  be  written  in  terras  of  z,  since  the  sliding 
velocity  is  constant. 


Then,  at  y  =  0 

d02  '  K2T  d*d! 
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where 


1/2 

I2  =  z  cos  uj^x  cos  a 2 z  dxdz 


(7) 


Pressure  Fluctuation  Due  to  Heat  Disturbance  -  Assume  that  the 
change  in  shape  (v)  due  to  thermal  deformation  can  be  approximated  as  a 
cosine  function  in  the  xy-plane: 
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v  =  C  cos  —  x  (8a) 

a 

and  in  the  yz-plane 

v  =  C  cos  —  z  (8b) 

b 

where  a  and  b  are  the  dimensions  of  the  contact  region  illustrated  in 
Figure  18  and  C  is  the  amplitude  of  the  perturbation.  This  change 
increases  the  deformation  in  the  center  of  the  contact  and  reduces  it 
near  the  edges. 


The  changes  in  radius  of  curvature  due  to  the  thermal,  deformation 
are  defined  as: 


I 


and 


d 


2 

v 


1 


in  body  l.  Similar  changes  in  radius  also  occur  in  body  2.  (The  t 
subscript  is  related  to  the  changes  in  geometry  due  to  thermal  distor¬ 
tion.)  Then  from  the  assumed  distorted  shape  (Equation  8a  and  8b),  the 
radii  can  be  related  to  the  original  Hertzian  contact  dimensions. 


(9a) 


(9b) 


The  same  result  can  be  obtained  if  the  thermal  deformation  is  expressed 
as  a  Fourier  series. 


The  average  change  in  curvature  [23]  in  the  central  region  of  the 
contact  (where  the  presssure  is  increased)  if  the  surfaces  were  free  to 
expand  due  to  the  heat  input  is 

92Vj  ( 1  +  v  i  )ce  i  0, 

~~2  +  —2  =  (T  n  ab/4 

3x  5z  I 

^V2  (i  +  V2)a2  .  °2 

T  2  +  a  2  =  IT  Tab/4 

4x  5  z  2 

These  equations  can  be  used  to  find  the  geometrical  constants  At  and  Bt. 
The  increased  force  (Pt)  necessary  to  keep  the  central  part  of  the 
surface  flat  can  be  found  from  Hertz  equation  [24]  noting  that  the 
positive  perturbation  width  is  half  the  original  contact  width  (a). 


(10a) 


(10b) 
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The  frictional  heat  generated  by  this  positive  change  in  force  at 
the  center  of  contact  is 


0  =  fP  V 

F  t 


(12) 


There  will  also  be  a  redaction  in  load  at  the  outer  edge  of  the  contact 
and  reduced  heat  generation  in  that  region,  but  that  will  not  lead  to  an 
unstable  condition. 


The  Thermoelastic  Stability  Condition  -  Substituting  equations  (5), 
(7),  and  (12)  into  the  heat  balance  equation  (1)  at  the  center  of  the 
contact,  an  expression  for  the  thermoelastic  stability  is  obtained. 
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The  above  equation  gives  the  relation  among  load ,  velocity,  coeffi¬ 
cient  of  friction  and  material  properties  for  the  onset  of  instability. 
If  the  materials  of  the  two  cylinders  are  the  same: 
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where  mt  and  r*t  are  geometric  coefficients  of  the  distorted  shape  from 
Hertz  analysis  [24].  Equation  (14)  gives  the  value  of  "a"  when  the  heat 
balance  at  the  center  of  the  contact  is  exactly  satisfied. 

This  can  be  translated  to  a  total  system  load  using  equation  (11) 
and  the  overall  geometry  of  the  contact  rather  than  the  perturbation 
geometry.  The  resulting  transition  load  (W)  between  stable  and  unstable 
operating  conditions  can  then  be  written  as: 


w 


43 


144  (A  +B) 

2  2 
E  (1-0 


(■n  mt^t  K  (  1  — v )  \  3 
fVma 


(15) 


Results  and 


Discussion  of  the  Technique  Used  for  the  Simplified  Case 


Equation  (15)  was  evaluated  for  the  case  of  a  crowned  steel  disk 
(body  1)  in  contact  with  a  cylindrical  disk  (body  2)  of  the  same 


*ial« 

,  The  specific  values  of 

the  parameters  are: 

K  = 

14.4  W/m*K 

Ri  - 

8.24 

mm 

v  = 

0.3 

Rl’  = 

53.2 

mm 

E  * 

206.8  x  109  N/m2 

R2  = 

31.7 

mm 

t  = 

17.3  x  106  /  °C 

r2’  = 

CO 

The  calculated  scuffing  loads  for  different  sliding  speeds  from  Equation 
(15)  are  plotted  in  Figure  19  with  a  range  of  values  for  friction 
coefficient  from  0.04.  to  0.2  illustrated.  A  lower  friction  coefficient 
(or  reduced  sliding  speed)  will  result  in  less  heat  generation  at  the 
interface  and  therefore  less  thermal  distortion.  As  a  result,  a  higher 
load  will  be  required  to  produce  enough  distortion  for  the  unstable 
heating  situation  to  exist;  that  is,  scuffing  failure  to  occur.  In  the 
same  way,  higher  values  of  friction  coefficient  (or  increased  speed)  will 
lead  to  lower  scuffing  loads. 


Also  shown  in  Figure  19  are  experimental  values  of  scuffing  loads 
taken  from  Durkee  [17].  These  data  are  for  two  pairs  of  steel  surfaces 
with  different  roughnesses:  "smooth"  pairs  which  consisted  of  a  0.025  urn 
cla  lower  surface  in  contact  with  a  0.05  um  cla  upper  surface,  and 
"rough"  pairs  which  consisted  of  a  0.25  um  cla  lower  surface  in  contact 
with  0.1  um  cla  upper  surface.  These  surfaces  were  lubricated  with 
naphthenic  mineral  oil.  The  experimental  scuffing  loads  were  strongly 
influenced  by  the  surface  roughness  of  the  mating  surfaces.  The  smooth 
pairs  displayed  scuffing  loads  approximately  one  order  of  magnitude 
larger  than  the  rough  pairs  at  the  same  sliding  speed.  Comparing  the 
experimental  results  to  the  corresponding  constant  friction  coefficient 
curves  from  the  analysis  indicates  a  value  of  f  less  than  0.1  for  the 
smooth  pairs  and  approximately  0.2  for  the  rough  surfaces. 

The  maximum  Hertzian  stress  for  the  elliptical  contact  is  shown  in 
Table  3  for  loads  from  200  to  1000  N  (45  to  225  lbs). 


Table  3  Contact  geometry  for  different  normal  loads 


Load , 

N 

Ha If -width 
of  contact, 

mm 

Hal f -length 
of  contact, 

mm 

Maximum 

Hertzian  Stress, 

GPa 

200 

0.548 

0.139 

1.26 

400 

0.691 

0.175 

1.58 

600 

0.791 

0.200 

1.81 

800 

0.870 

0.220 

2.00 

1,000 

0.938 

0.237 

2.15 

The  elastohydrodynamic  film  thickness  for  these  conditions  (high  contact 
stress  and  pure  sliding)  will  be  extremely  thin  and  thus  contact  of  the 
surface  asperities  will  most  certainly  occur.  This  would  explain  the 
higher  apparent  friction  coefficient  for  the  rough  surfaces.  An  increase 
in  roughness  would  be  accompanied  by  more  asperity  contacts,  smaller 
percentages  of  the  load  supported  by  the  fluid  film,  and  higher  friction 
force.  EHD  film  would  also  explain  the  apparent  reduction  in  friction 
coefficient  for  the  smooth  surfaces  as  the  speed  is  increased.  Higher 
speeds  will  generate  a  thicker  film  and  a  net  decrease  in  friction  force. 
Therefore,  to  predict  the  boundary  of  the  stable  region,  that  is  the 
scuffing  load,  an  independent  expression  for  friction  coefficient  must  oe 
used. 


An  estimate  of  the  friction  coefficient  [25]  in  the  mixed 
lubricated  regime  is: 


f  =  0.0127  log10(- 


3.17  x  108  W/L 
u ( V j -  V2)(Vj+  V2)‘ 


(16) 


where  Vj  and  V2  are  in  in/s,  W/L  load  per  unit  width  in  lb/in,  and  u  the 
absolute  viscosity  in  centipoise.  In  the  elliptic  contact  of  smooth 
surfaces,  an  equivalent  load  W/L  is  used  in  this  equation.  The 
equivalent  load  per  unit  length  is  that  load  which  produces  the  same 
Hertzian  contact  stress  in  the  line  contact  of  two  cylinders  as  that  in 
the  elliptic  contact  case. 

In  the  case  of  rough  lubricated  surfaces,  the  coefficient  of 
friction  at  V  =  4.15m/s  is  chosen  as  0.14  from  experimental  results  [17]. 
The  coefficient  of  friction  at  other  speeds  is  calculated  in  proportion 
to  f  =  0.14  by  means  of  Equation  (16)  with  constant  values  of  W/L  and  11. 
The  coefficient  of  friction  thus  obtained  is  shown  in  Table  4. 
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Table  4  Friction  coefficient  [21]  for  sliding  surfaces 


Speed 

2  m/s 

3  m/s 

4  m/s 

5  m/s 

6  m/s 

7  m/s 

f  (smooth) 

0.068 

0.060 

0.055 

0.051 

0.048 

0.045 

f  (rough) 

0.  158 

0. 148 

0. 140 

0.134 

0.129 

0.  126 

Figure  20  illustrates  the  agreement  between  the  measured  scuffing 
loads  and  the  theoretical  predictions  of  thermally  unstable  contact 
regions.  In  general  the  agreement  is  quite  good;  in  most  cases,  the 
differences  between  calculated  and  measured  limiting  speeds  are  below 
30%.  The  difference  at  high  load  is  large,  reaching  72%  at  940N  load. 
The  cause  could  be  that  the  friction  coefficient  in  that  case  (0.068)  may 
be  too  low.  Lubrication  is  poor  at  low  speed  and  high  load  and  Equation 
(16)  may  not  predict  the  value  under  these  conditions. 


Experimental  Testing  of  Failure  Due  to  Scuffing 
Apparatus  Design 

The  main  parts  of  the  experimental  equipment  are  two  disks  in 
contact.  The  upper  disk  is  45  mm  in  diameter  with  a  crown  radius  of 
140  mm.  The  lower  disk  is  cylindrical:  60  mm  in  diameter  and  85  mm  in 
length.  This  cylinder  is  supported  by  two  ball  bearing,  and  is  driven  by 
a  3  HP  variable  speed  motor  through  a  V-belt  drive.  The  runout  of  the 
cylinder  is  below  13  um.  A  belt  drive  is  used  to  get  a  maximum  speed  of 
3000  RPM.  The  highest  speed  used  in  the  experiment  has  been  2100  RPM.  A 
reduction  belt  drive  is  used  at  speeds  below  1000  RPM  to  get  sufficient 
torque  by  changing  the  pulleys  of  the  belt  drive.  The  upper  disk  is  also 
supported  in  ball  bearings  located  in  the  case  with  a  load  arm  on  either 
end.  In  the  present  experiments,  the  upper  disk  is  held  stationary  by 
the  friction  measurement  device.  Modifications  may  be  made  utilizing 
measuring  and  controlling  devices,  to  allow,  the  upper  disk  to  rotate 
during  the  tests. 

After  each  test,  both  the  surfaces  of  the  upper  and  lower  disks  are 
either  damaged  or  worn.  New  regions  on  the  two  disks  are  brought  into 
contact  for  the  next  test  by  moving  the  upper  disk  axially  a  certain 
distance  along  the  cylinder.  Because  the  upper  disk  does  not  rotate,  it 
can  be  moved  through  a  small  angle  to  get  a  new  part  of  surface  on  the 
upper  disk  into  the  contact  region. 

The  load  is  applied  by  hanging  weights  on  the  end  of  the  load  arm. 
A  counterbalance  load  is  fixed  on  the  other  side  of  the  load  arm. 
Weights  are  suspended  on  this  end  of  the  arm  so  as  to  counterbalance  the 
weight  of  the  upper  disk  and  the  load  arm  on  the  other  side.  Thus,  there 
is  no  force  between  the  two  disks  when  no  weight  is  hung  on  the  load  arm. 
The  maximum  force  applied  on  the  contact  is  2100  N  (465  lb.)  when  weights 
of  147  N  (33  lb.)  are  hung  on  the  load  arm. 
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The  rotation  speed  is  measured  by  means  of  a  toothed  gear  and  a 
transducer.  The  gear  with  60  teeth  is  fixed  on  one  end  of  the  shaft  of 
the  lower  cylinder.  The  displacement  transducer  is  placed  near  the  gear 
and  the  output  is  sent  to  a  frequency  counter  which  indicates  the 
rotation  speed. 

Solenoids  are  placed  under  the  load  arm  to  facilitate  releasing  the 
load  by  pushing  the  j.oad  arm  upward  at  the  instant  when  surface  failure 
occurs.  The  solenoids  can  be  triggered  by  the  temperature  rise,  or  the 
friction  increase,  or  the  vibration  produced.  The  vibration  level  as 
measured  via  an  accelerometer  was  found  to  be  the  best  indicator  of 
scuffing  failure. 

Temperature  in  the  contact  area  of  the  two  disks  during  sliding  is 
sensed  by  an  infrared  thermo-monitor  (Vanzetti  model  3025)  and  is 
displayed  on  an  oscilloscope.  The  response  of  the  infrared  thermometer 
is  fast  and  no  temperature  gradient  is  produced  as  in  the  ca • e  of 
thermocouple,  so  it  is  more  effective  in  monitoring  transient  tempera¬ 
ture.  The  spot  size  of  the  focused  light  is  about  1.-0  mm  and  the 
distance  of  focus  is  38  mm.  Unfortunately,  the  infrared  thermometer  can 
only  sense  relatively  high  temperature;  the  minimum  temperature  sensed  is 
about  80°  C.  The  purpose  of  this  research  is  to  investigate  the 
temperature  rise  during  scuffing  of  a  sliding  contact;  this  temperature 
is  typically  above  80°  C.  If  the  contact  temperature  is  below  80°  C,  the 
scuffing  regime  has  not  been  reached. 

The  friction  force  produced  in  the  contact  during  sliding  is 
measured  by  the  device  which  is  fixed  on  the  shaft  of  the  upper  disk. 
The  friction  force  tends  to  rotate  the  upper  disk  and  with  it  a  bar  which 
is  clamped  on  the  shaft  of  the  upper  disk.  The  bar  compresses  a  spring 
and  changes  the  small  gap  between  the  surfaces  of  an  eddy  current 
proximity  sensor  and  the  bar,  and  an  output  proportional  to  the  traction 
force  is  detected. 

The  bolt  on  the  side  of  the  spring  is  used  to  adjust  the  gap  between 
the  transducer  and  the  bar  and  the  preload  of  the  spring.  The  diameter 
of  the  bolt  hole  is  made  bigger  than  that  of  the  bolt,  so  the  bar  is  free 
to  rotate  without  constraint.  The  hook  on  the  bar  is  used  to  attach  a 
wire  for  hanging  weights  during  calibration.  The  wire  passes  through 
pulleys  which  are  clamped  on  the  base  of  the  test  rig  only  at  the  time  of 
calibration.  On  the  other  end  of  the  wire,  weights  of  known  magnitude 
are  hung.  At  the  same  time  readings  of  voltage  on  the  oscilloscope  ate 
taken.  Thus,  the  relation  between  force  and  voltage  reading  can  be 
found . 


The  cylinder  and  the  disk  are  lubricated  by  a  wick  system.  A  basin 
is  placed  under  the  lower  cylinder  with  oil  poured  in  it.  Some  pieces  of 
felt  are  fixed  on  the  bottom  of  the  basin.  The  felt  touches  the  lower 
surface  of  the  cylinder.  When  the  cylinder  turns,  sufficient  oil  is 
carried  on  the  surface  of  the  cylinder  into  the  contact  region. 

Test  Procedure  -  During  the  test,  the  speed  of  the  lower  cylinder  is 
adjusted  to  the  desired  rotating  speed.  A  certain  amount  of  load  is 
added  to  the  load  arm  and  the  test  is  run  for  this  load  from  the 


beginning  Co  the  end.  If  scuffing  failure  occurs,  temperature,  friction 
and  vibration  increase  steeply  and  simultaneously.  At  this  point  the 
solenoids  are  activated  and  the  load  arm  is  pushed  up  to  separate  the  two 
contact  surfaces  and  the  test  stops.  If  failure  does  not  occur,  the  test 
is  run  for  about  60  seconds  and  then  is  stopped.  A  sequence  of  tests  are 
run  in  which  the  load  is  either  less  than  the  previous  test,  or  larger, 
depending  on  whether  there  was  failure  or  not.  The  tests  are  continued 
until  the  threshold  of  scuffing  load  is  obtained. 

The  scuffing  tests  were  run  in  groups  according  to  the  combination 
of  lubricant  and  disk  material  used.  Each  group  consisted  of  approxi¬ 
mately  35  scuffing  tests  which  were  run  over  the  load  and  speed  ranges  of 
the  apparatus.  The  following  is  a  listing  of  the  test  groupings: 

Test  Group  1 

'»  lubricant:  XRM-109  synthetic  paraffinic  oil,  viscosity  =  447 
cStokes  at  40  C 

**  roller  material:  52100  alloy  steel,  hardness  Rc  60 

*  calculated  specific  film  thickness  (X  ratio):  5  to  8 
Test  Group  2 

*  lubricant:  Mobil  SHF-61,  synthetic  base  stock,  viscosity  =  28 
cStokes  at  40  C 

**  roller  material:  52100  alloy  steel,  hardness  Rc  60 

calculated  specific  film  thickness  (X  ratio):  0.4  to  1.3 
Test  Group  3 

*  lubricant:  Mobil  SHF-401,  synthetic  base  stock,  viscosity  =  443 
cStokes  at  40  C 

3  roller  material:  52100  alloy  steel,  hardness  Rc  60 

n  calculated  specific  film  thickness  (X  ratio):  4  to  12 
Test  Group  4 

lubricant:  Mobil  SHF-401,  synthetic  base  stock,  viscosity  =  443 
cStokes  at  40  C 

**  roller  material:  440C0  stainless  steel,  hardness  Rc  60 

*>  calculated  specific  film  thickness  (X  ratio): 


7  to  20 
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The  following  parameters  were  either  determined  or  set  prior  to  a 
test:  geometry  of  contact,  surface  roughness  of  contact,  load,  speed, 

lubricant,  and  roller  material.  The  parameters  measured  during  a  test 
were:  frictional  force,  vibration,  temperature,  and  time  to  failure. 

The  sequence  of  tests  within  a  test  group  generally  depended  upon 
whether  failures  had  occurred  during  the  first  several  tests.  If  a 
failure  had  occurred,  either  load  or  speed  was  reduced  in  subsequent 
tests  until  a  non-failure  was  observed.  If  a  failure  did  not  occur 
initially,  load  or  speed  was  increased  in  each  subsequent  test  until 
failure  occurred.  Each  test  was  run  on  a  new,  un-worn  spot  on  the  fixed 
disk  and  on  a  new  un-worn  Crack  on  the  rotating  disk.  In  this  way,  a 
failure  boundary  was  defined.  The  choice  as  to  whether  to  vary  load  or 
speed  depended  upon  the  region  of  the  theoretical  fV  curve  being 
explored.  At  high  speeds  and  low  loads,  the  curve  predicts  relative 
insensitivity  to  further  changes  in  speed.  At  high  loads  and  low  speeds, 
the  curve  predicts  insensitivity  to  further  changes  in  load. 

Experimental  Results 


The  theoretical  scuff  boundary  curves  presented  herein  were  used  to 
evaluate  the  results  of  the  experiments  and  to  correlate  them  with  the 
theory.  The  theory  would  be  shown  to  be  ideally  valid  if  all  non- 
failures  occurred  in  both  plots  on  the  side  of  the  boundary  closest  to 
the  ordinate.  Conversely,  scuffing  failures  would  ideally  occur  on  the 
opposite  side  of  the  boundary. 

The  intent  of  each  test  group  was  to  evaluate  the  response  of  the 
model  for  a  different  lubricant-metal  combination;  therefore,  the  results 
are  presented  by  test  group.  A  plot  is  given  for  each  test  group:  load 
versus  heat  generation,  0  =  fVL.  In  each  of  the  figures,  failures  are 
indicated  by  circles  and  nonfailures  are  indicated  by  squares.  The 
theoretical  boundary  is  shown  by  a  solid  line. 

The  data  of  group  1  (Figure  21)  correlates  with  the  theoretical 
prediction  of  the  failure  boundary,  in  that  most  of  the  experimental  test 
results  occurred  on  the  side  of  the  boundary  for  which  they  were 
predicted  to  occur.  The  uncertainty  band  of  t  10%  is  included  in  the 
f igure. 

Group  2  results  (Figure  22)  show  a  marked  difference  from  the 
results  of  Group  1  in  that  no  correlation  with  the  predicted  failure 
boundary  is  apparent.  The  difference  between  the  two  test  groups  was  the 
contact  lubricant  used.  SHF-61  is  much  less  viscous  than  XRM-109.  Re¬ 
duction  of  lubricant  viscosity  reduces  the  film  thickness  in  the  contact 
and  alters  the  coefficient  of  friction,  resulting  in  scuffing  failure 
under  less  severe  loads  and  velocities.  For  this  group  of  tests  X  lies 
in  the  range  of  0.4  to  1.3.  The  results  indicate  that  the  theory  does 
not  adequately  predict  the  scuffing  boundary  in  this  range  of  calculated 
specific  film  thickness.  In  addition,  the  calculated  X  ratios  indicate  a 
high  degree  of  asperity  contact,  suggesting  that  the  heat  generation  in 
these  tests  is  due  to  soldi  contact,  rather  than  shear  in  the  lubricant. 
The  theoretical  model  is  not  affected  by  this  distinction,  however,  be¬ 
cause  it  is  formulated  on  the  basis  of  an  overall  friction  coefficient  f. 
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Since  the  data  of  Group  2  do  not  correlate  with  the  prediced  failure 
boundary,  the  question  arises,  can  a  failure  boundary  be  determined  for 
this  set  of  data?  One  may  speculate  that  the  scuffing  boundary  curve  of 
the  model  could  be  modified  empirically  to  account  for  a  change  in  the 
calculated  specific  film  thickness.  It  is  apparent  that  the  actual 
failure  boundary  is  less  than  the  predicted  boundary.  This  implies  a 
modification  of  the  coefficient  of  friction  in  the  model  formulation  to 
suitably  account  for  lower  calculated  specific  film  thickness.  It  is 
also  possible  that  in  this  lower  specific  film  thickness  regime  another 
failure  mechanism  is  present  —  namely  that  described  by  Snidle,  Rossides 
and  Dyson  [26]. 

Test  results  for  Group  3  are  similar  to  those  of  Group  1,  as 
expected,  since  the  lubricants  are  very  similar.  Figure  23  illustrates 
that  the  results  correlate  with  the  theoretical  failure  prediction. 

The  objective  of  the  experiments  of  Group  4  (Figure  24)  was  to  test 
the  response  of  the  model  to  change  in  thermal  properties  (K,a),  although 
the  material  properties  (E,\>)  also  vary  when  the  disk  material  is 
different.  The  tests  were  run  in  the  high  calculated  X  ratio  regime, 
since  correlation  with  theory  had  already  been  observed. 

The  results  show  correlation  with  the  theoretical  failure  boundary, 
in  support  of  the  theoretical  model.  A  comparison  between  results  of 
Group  3  and  4,  in  which  the  disk  material  was  the  only  significant 
variable,  shows  that  the  experimental  results  do  not  answer  the  quesiton 
as  to  whether  there  is  a  difference  or  not. 

There  exist  two  possible  explanations  for  this  result.  One 
explanation  lies  in  the  possibility  that  the  friction  force  measurement 
was  not  sensitive  enough  to  detect  a  difference.  An  interesting 
experimental  observation  clouds  this  conclusion,  however.  As  described 
earlier,  a  sequence  of  tests  was  run  at  either  constant  load  or  speed  to 
define  the  experimental  failure  boundary.  This  is  observable  in  each  of 
the  plotted  figures  if  one  follows  a  line  of  constant  load,  for  example, 
across  a  plot.  In  many  cases,  despite  an  incremental  increase  in  speed, 
a  jump  in  the  magnitude  of  0  occurred  across  the  failure  boundary.  The 
jump  was  large  compared  to  the  instrument  resolution  capability  and  it 
indicates  that  an  increase  in  sensitivity  might  not  increase  the 
definition  of  the  failure  boundary,  and  that  the  scuffing  boundary  might 
be  thought  of  as  a  band  rather  than  a  line,  as  predicted  by  theory. 

The  second  explanation  is  that  the  response  of  the  model  to 
variation  in  thermal  and  material  properties,  although  theoretically 
predicted,  may  be  masked  in  actual  observation  by  other  mechanisms  or 
processes  not  adequately  addressed  by  the  model. 

Time  Dependency  of  Failure 

The  time  to  failure  does  not  appear  to  correlate  to  any  particular 
quantity  over  the  test  groups,  but  is  a  trend  observable  within  a  test 
set.  That  is,  within  a  set  of  tests  in  which  the  paired  surface  rough¬ 
nesses  did  not  change,  a  failure  generally  occurred  after  longer  period 
of  running  (or  not  at  all)  if  the  load  and/or  speed  was  reduced.  The 
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Results  of  Test  Group  1:  lubricant  XRM-109,  disk  material 
52100  steel,  5  <  X  <  8;  circles  indicate  scuffing  failures 
squares  indicate  non-failures. 
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Figure  22.  Results  of  Test  Group  2:  lubricant  SHF-61,  disk  material 
52100  steel,  0.4  <  X  <  1.3;  circles  indicate  scuffing 
failures,  squares  indicate  non-failures. 
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Figure  23.  Results  of  Test  Group  3:  lubricant  SHF-401,  disk  material 
''  52100  steel,  4  <  X  <  12;  circles  indicate  scuffing  failures, 

squares  indicate  non-failures. 


I 

I 


CONTACT  LOAD.  KG 


55 


Figure  24,  Results  of  Test  Group  4:  lubricant  SKF-401,  disk  material 
440C  stainless  steel,  7  <  X  <  20;  circles  indicate  scuffing 
failures,  squares  indicate  non-failures. 
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actual  test  duration  and  load/speed  relationship  was  not  consistent 
between  test  sets. 

This  general  trend  was  observed  in  test  Groups  1,3  and  4.  Test 
duration  in  Group  2,  in  which  the  calculated  X  ratio  was  significantly 
lower,  did  not  follow  the  described  trend.  Failures  generally  occurred 
within  the  first  minute,  or  not  at  all.  Most  failures  occurred  in  the 
first  few  seconds  of  running. 

These  observations  are  in  agreement  with  Durkee's  research  [18],  in 
which  scuffing  experiments  were  run  for  two  ranges  of  surface  finish. 
Durkee  observed  a  time  dependency  of  failure  for  the  "smoother"  tests  and 
immediate  failure  for  the  "rough"  tests.  Although  Durkee  does  not 
present  specific  film  thickness  values,  it  may  be  inferred  that  the  tests 
run  on  smooth  surfaces  correspond  to  larger  X  ratios,  relative  to  the 
tests  run  on  rough  surfaces. 

Wear  Observations 


Throughout  the  scuffing  experiments,  wear  was  visible  on  both 
cylinders  for  tests  in  which  failure  did  not  occur.  This  was  an 
unexpected  result,  since  calculated  specific  film  thickness  predicted 
little  or  no  asperity  contact  in  test  Groups  1,2  and  4.  The  disparity 
between  observed  wear  and  X  ratio  calculations  is  attributed  to  the 
calculation  model,  which  does  not  consider  heating  due  to  slip  in  the 
contact . 

The  weld  patterns  left  on  the  lower  cylinder  from  scuffing  failure 
were  observed  to  form  a  track  which  was  slightly  helical  about  the 
cylinder  axis,  suggesting  that  failure  occurred  at  one  point  in  the 
contact  region,  and  traversed  the  contact  as  material  was  lost  to  the 
rotating  cylinder.  This  correlates  to  the  moving  hot  spots  described  in 
studies  by  Dow  13.  The  dimensions  of  the  wear  markings  left  on  the  upper 
cylinder  when  a  failure  did  not  occur  were  measured  and  found  to  compare 
closely  with  the  calculated  Hertzian  dimensions. 

Effects  of  Exprimentai.  uncertainty 

Although  the  frictional  force  measurement  contained  a  relatively 
large  uncertainty,  the  effect  of  this  is  minimized  in  light  of  the  most 
important  determination  —  that  is,  prediction  of  failure.  How  near  or 
far  away  a  test  point  lies  to  the  predicted  boundary  is  of  secondary 
importance.  In  most  cases,  an  uncertainty  of  ±10%  in  the  frictional 
force  does  not  alter  which  side  of  the  boundary  a  test  point  lies.  It  is 
interesting  to  not  that  relatively  few  data  points  were  actually  found  to 
lie  in  proximity  to  the  boundary,  although  considerable  effort  was  made 
to  obtain  data  in  the  boundary  region. 
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